Filtering of matter wave vibrational states via spatial adiabatic passage 
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We discuss the filtering of the vibrational states of a cold atom in an optical trap, by chaining 
this trap with two empty ones and controlling adiabatically the tunneling. Matter wave filtering is 
performed by selectively transferring the population of the highest populated vibrational state to the 
most distant trap while the population of the rest of the states remains in the initial trap. Analytical 
conditions for two-state filtering are derived and then applied to an arbitrary number of populated 
bound states. Realistic numerical simulations close to state-of-the-art experimental arrangements 
are performed by modeling the triple well with time dependent Poschl-Teller potentials. In addition 
to filtering of vibrational states, we discuss applications for quantum tomography of the initial 
population distribution and engineering of atomic Fock states that, eventually, could be used for 
tunneling assisted evaporative cooling. 
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I. INTRODUCTION 

Ultracold atoms trapped in optical potentials [l-Q, 
e.g., microtrap arrays or optical lattices, have attracted 
considerable attention since they fulfill all the basic re- 
quirements for quantum information processing [J]. In 
fact, neutral atoms in dipole trap arrays with short- 
range interactions such as s-wave scattering of bosons [5( 
or state selective long range interactions such as dipole- 
dipole interactions |6J], do not experience intrinsic limi- 
tations in their scalability. Thus, a quantum register of 
about a hundred qubits has been reported recently in 
a two-dimensional (2D) optical microtrap array [3[. 2D 
optical microtrap arrays present two characteristic fea- 
tures: (i) the simplicity to achieve single-site addressing 
since the trap separation distances can range from single 
/xm up to ~ 100 /xm, in any case being larger than the 
resolution limit; and (ii) the freedom to move indepen- 
dently sets of traps or even individual traps to control 
the tunnelling. However, cooling a single atom down to 
the lowest vibrational state of an optical microtrap array 
is still a challenging issue and, for some physical real- 
izations, represents one of main experimental limitations 
to perform quantum computations with optical micro- 
traps. Thus, developping novel techniques to determine 
and, eventually, to engineer the population distribution 
in optical microtraps is a focus of present research [3| . In 
this context, we here propose to make use of the spatial 
adiabatic passage technique Q to achieve this goal. 

Spatial adiabatic passage consists in adiabatically fol- 
lowing a spatial dark state whose spatial profile is de- 
termined by the tunneling interaction between neighbor- 
ing traps [7| and, in fact, is the matter wave analog of 
the well known quantum optical Stimulated Raman Adi- 
abatic Passage (STIRAP) technique Q. Here, we will 



take profit of the fact that tunneling rates between traps 
strongly depend on the vibrational state under consider- 
ation to perform state selective adiabatic passage leading 
to (i) filtering of vibration states, (ii) quantum tomogra- 
phy of the initial population distribution, and (iii) engi- 
neering of atomic Fock states that, eventually, could be 
used for tunneling assisted evaporative cooling. 

The article is organized as follows. In Section II we 
introduce the physical model consisting in a single atom 
in three identical Poschl-Teller type potentials Q with 
time varying position of the trap centers. The filtering 
protocol is presented in Section III and analytical con- 
ditions for two-state filtering are derived and compared 
to numerical simulations. Section IV focuses on the ap- 
plication of the filtering protocol to multiple-states for 
quantum tomography and quantum engineering of Fock 
states. Finally, the conclusions are summarized in Sec- 
tion V. 



II. MODEL 

We study the dynamics of a single cold neutral atom 
of mass to in a ID triple well potential, see Figs.QJa) and 
(b), described by the Schrddinger equation: 






1m dx 2 



i>(x,t), (1) 



where the wells are modeled by three identical Poschl- 
Teller (PT) type potentials |j|: 



V{x,t) = J2 Voisech 2 \y/2(x-Xi(t))/ 



i=L,M,R 



(2) 



Voi is the potential depth for the i-th trap, and Xi (t) 
defines the position of its center at time t. 2a is the 



width of the PT potential and tu x = h/ma 2 is the trap- 
ping frequency. Note that we assume here PT potentials 
since they can be used to model very accurately Gaus- 
sian potentials (see Appendix A), obtained with dipole 
traps built-up by focusing a laser beam. In this case the 
parameter a corresponds to the waist wo of the Gaussian 
light beam. In addition, PT potentials provide analytical 
expressions for their energy eigenvalues and eigenstates, 
see Sec. Ill B. Throughout the paper, we will use dimen- 
sionless units for time, tw x , space, x/a, and potential am- 
plitude, Vo/huj x . Initially, at t inil we will assume that the 
neutral atom is distributed among the vibrational states 
of the left trap, while the other two traps are empty. 



III. TWO-STATE FILTERING 

A. Basic idea 

Our proposal for the filtering of vibrational states is 
based on an adiabatic transport process |7[ between the 
two extreme traps that resembles the well known quan- 
tum optical Stimulated Raman Adiabatic Passage (STI- 
RAP) technique [8J. Although we will accurately inves- 
tigate the filtering of vibrational states by numerically 
solving the Schrddinger Eq. (1), for simplicity, in the fol- 
lowing lines we will illustrate the basics of our proposal by 
restricting the dynamics into the Hilbcrt space spanned 
by the ground and the first excited vibrational states of 
each trap. In this finite basis, the Hamiltonian of the 
system can be roughly approximated to H = Ho Hi 
with: 



o n% M {t) o 
H n = h\ n^(t) o n™ R (t) 
o n^ R (t) o 



(3) 



where f2Jf is the tunneling rate between two adjacent 
traps i and j with i,j — L,M,R denoting left, middle 
and right, respectively and n = 0, 1 refers to the ground 
or the first excited vibrational state, respectively. Note 
that we have split Hamiltonian H into the direct sum of 
Hq and Hi and, therefore, we have assumed that the en- 
ergy separation between the ground and first vibrational 
state of each trap is large enough to avoid crossed tunnel- 
ing between ground and excited states of different traps. 
This approximation could fail for vibrational states close 
to the continuum, such as the vibrational analogues of 
Rydberg states, where the energy spacing between differ- 
ent levels becomes relatively small. The latter scenario 
is out of the scope of this paper. 

After diagonalization of the two Hamiltonians given 
in Eq. (3), one ends up, in particular, with two energy 
eigenstates that only involve vibrational states of the two 
extreme traps: 

\D n {9 n )) = cos6 n \n) L - sm9 n \n) R withn = 0, 1(4) 

where the mixing angle, 6 n , is defined as tan#„ = 
Q^ M /Q,^f R , States \D n (9 n )) are known as spatial dark 
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FIG. 1: (Color online) (a) Temporal variation of the po- 
sitions of the PT traps centers, where do is the minimum 
trap separation and T is the time delay between the two ap- 
proaching sequences, (b) Illustration of the spatial profile 
of the PT potentials at the three different times correspond- 
ing to the vertical lines in (a). Each isolated PT trap has 
depth of Vo — —20fkd x and supports only 4 vibrational states, 
(c) Temporal variation of the tunneling rate fl^ M (Q^f R ) be- 
tween left and middle (middle and right) traps for the ground 
(n = 0) and first (n = 1) excited states. 6\ is the mixing 
angle for the first excited states. Temporal variation of the 
population of the ground and first excited states of the left (d) 
and right (e) traps assuming the following initial distribution 
Po(t = 0) = Pi {i = 0) = 1/2, respectively. 



states 0- State-selective adiabatic passage of matter 
waves between the two extreme traps will consist in adia- 
batically following one of the two energy eigenstates (4) , 
typically the one with highest energy, by the smooth vari- 
ation of the tunneling rates, while for the other one the 
transport process is inhibited. 



B. Poschl— Teller tunneling rates 



Therefore, 



For a single PT potential © of depth V = -s (s + 1) 
(s > 0) there is an analytical solution for their energy 

cigenstates in terms of the associated Lcgendrc PJ n_s) 
polynomials as follows: 



<t> n , s (x) = iV„, s PJ"- s > [tanh(V2 (x - x t ) /a) 
with normalization constant: 

2 l/4 



N n . = 



,1/2 



(s - n) 



T(2s-n + l) 
r(n + l) 



(5) 



(G) 



where s = y/\V \ + 1/4— 1/2, n = 0, 1, . . . , N max numer- 
ates the bound states with N max being the integer part 
of s, and E n , s /hu> x — — (s — n) gives the energy of the 
n-th state. T is the Gamma function. 

The tunneling rates fl n . s between two identical PT 
potentials can be found by determining the energy dif- 
ference between the symmetric 0+ s and antisymmet- 
ric 0~ s eigenstates namely il n . s = (0+ 8 |if| 0+ s ) - 
{K,s\H\ct>- >s ) = E+ s - E- s . We take </>±' s = 
(4> % n a ±9^ s ) /\/2 where <j$ a corresponds, see Eq. (5), 
to the localized state <f> ntS of either the trap i or its neigh- 
bour j. For n = 0, the Gram-Smith orthonormalization 
procedure (see Appendix [B| provides very accurate an- 
alytical expressions for the ground state tunneling rate, 
f^o.s- For the exited states, only approximate analytical 
solutions are possible. The Holstein-Herring method, see 
[Hj[ll|], yields for two identical traps: 



n n>a (d) = e+s- = 



-V «.(*))' 



1-2J? (&,&))'& 

(7) 
where d is the distance between the trap centers. 

At large distances between the traps, the denominator 
in ([7J) rapidly approaches unity and the main behavior of 
the tunneling rate J7„ jS (d) is given predominantly by the 
numerator, for which by applying the recurrence relations 
for associated Legcndre functions one could write 



Q n , 8 (d) ~ -V(CW) 



x=d/2 

2cf>i l (x)N n , s [V(l-i 2 )i J i"- s)+1 + 



+ (n-s)xP^ n - s) 



x— tanh d/2 



(8) 



For tanh d/2 — > 1 the above expression could be further 
simplified 

0„, 8 (d) - 2 (<^ s (d/2)) 2 (s - n) = S„, s e-("- n ) d , (9) 
where 



B n . s — 



r(2s-n+l) f2s(s + l) 



r(n+l) \T(s-n) 



ftn.s (d) 
(d) 



n 



An 



n—l,s 



with 



A 



(s-n) 2 



(2s- n+ l)n 



(11) 



(12) 



From (9), it is clearly shown that for a fixed poten- 
tial depth (fixed s) the tunneling rate £l n ,s (d) increases 
with the energy level n and decreases with the distance 
d. In contrast, the tunneling rate ratio for two consecu- 
tive levels, Eq. (11), exhibits inverse dependence, i.e. it 
decreases with the energy level n and increases exponen- 
tially with the distance d. 



C. Conditions for two-state filtering 

For the transfer process we will assume that the posi- 
tion of the trap center for each of the traps can be varied 
at will to temporally control the tunneling interaction 
|3,Q. In this case, the adiabatic transport will consist in 
approaching and separating the traps in a counterintu- 
itive sequence, see Fig.flja), with typical spatial profiles 
given in Fig. Hfb). The two empty traps, right and mid- 
dle ones, are approached and separated first and, with 
an appropriate time delay T, left and middle traps are 
approached and separated. The motional sequence of the 
traps, see Fig. Q3a), * s engineered in such a way that the 
time variation of the tunneling rates between two ad- 
jacent traps resembles a Gaussian profile, see Fig. QJc). 
With this aim and taking into account the explicit de- 
pendence of the tunneling rates with the distances, see 
Eq. (9), we fix xu = and take the following temporal 
variation for the outermost trap positions 



(x L - x M ) 
a 

(x R - x M ) 




t-l\ -2L 



\auj x 



a 



T 

f+ 2 



_vo_ 
aui a 



^ ,l;n 

a 



with do being the minimum separation distance between 
the outermost (either left or right) and the middle trap 
achieved at time t = ±T/2, respectively, vq gives the 
modulus of the velocity of the outermost traps at large 
separation distances. 

Since we are considering identical traps, the tunnel- 
ing couplings between right and middle and between left 
and middle traps will follow the same dependence with 
the trap distance, fl^ M (d) = ^ R {d) = Cljd). In this 
case, the "global" adiabaticity condition 0, @ for spatial 
adiabatic passage reads: 



(10) ^(^(do)) 2 + (^ n IR (d )) 2 T = n n (d )T > 10 (14) 



where T is the characteristic time for the adiabatic pas- 
sage process. 

Since we are interested in state-selective atom trans- 
fer, the goal of the filtering protocol will be that the 
atomic population initially in the highest vibrational 
state of the left trap follows adiabatically the spatial dark 
state, |-Di(#i)) ending in the right trap, while the popu- 
lation initially distributed in the lower levels of the left 
trap remains there during the whole process. Therefore, 
for the two-state filtering case, the counterintuitive mo- 
tional sequence of the traps should be performed fulfill- 
ing: fli(d )T > 10 and fl (d )T < 10. Moreover, the 
filtering protocol requires also to inhibit the direct trans- 
fer between neighbouring traps of the population initially 
in the ground state of the left trap, i.e., H,o(do)T <C 1. 
Therefore, the necessary condition for two-state filtering 
reads: 



»i(do) 



>10. 



(15) 



Note also that in order to avoid the direct coupling 
between the outermost traps, it is also required that 
fl n (2do)T <C 1 which implies: 



»i(do) 
fii(2do) 



>10. 



(16) 



For the filtering sequence shown in Fig. Q] (a)-(c) with 
6\ varying from to 7r/2, we have choosen parame- 
ters such that both conditions (fT5j) and (fTB"|) are fulfilled 
and, therefore, one expects that the filtering protocol 
succeeds. Figs. QJd) ancl IH e ) plot the temporal varia- 
tion of the population distribution of the left and the 
right traps, respectively, by integrating the correspond- 
ing Schrodingcr Eq. (1) with the initial population dis- 
tribution P L (t = 0) = Pf(t = 0) = 1/2. At the end of 
the process, Pq = P^ = 1/2, which confirms the validity 
of the filtering protocol. 

In the following, we will investigate the robustness 
of the filtering protocol under variations of the param- 
eters. With this aim, we plot in Fig. [21(a)- (c) curves 
Qi(d )T = 10 (solid blue) and fl (d Q )T = 1 (dashed 
red) in the parameter plane (do /a, uj x T) for three differ- 
ent values of the potential depth s = 2, 3 and 4. The 
dotted green curve corresponds to Qi(2d )T = 1. The 
grey region defines the parameter domain for which both 
conditions (15) and (16) are fulfilled. Note from Fig. [2] 
(a)-(c) that even for small values of s, the parameter do- 
main where the filtering protocol should succeed is lim- 
ited by condition (15). To confirm the previous predic- 
tions, we have performed numerical simulations of the fil- 
tering protocol integrating the Schrodingcr equation for 
s — 4. Fig. EJd) shows the contour plot of the fidelity at 
the end of the filtering process, defined as F = P L + P^ 
for an initial population distribution of Pq = Pf = 1/2. 
For this case, the previously derived filtering conditions 
(fTo]) and (fTB")) assure that the fidelity of the process is 
above 0.99. We have numerically checked the validity of 
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FIG. 2: (Color online) Curves fii(d )T = 10 (solid blue) and 
Qo(do)T = 1 (dashed red) in the parameter plane (do /a, u x T) 
for three different values of the potential depth (a) s — 2, (b) 
s = 3 and (c) s = 4. In all three cases, the dotted green 
curve corresponds to Qi(2do)T — 1. The grey region defines 
the parameter domain for which the filtering conditions (|15[) 
and (|16|) are fulfilled, (d) Contour plot of the fidelity (see 
text) of the filtering process for s — 4 obtained numerically 
by integration of the Schrodinger equation with the temporal 
variation of the traps centers given by (|13[) . The cross in (d) 
marks the parameter setting used in Fig. [TJ 



the derived filtering conditions for a wide set of parame- 
ters. 



D. Extension to two excited states 

The generalization of the filtering protocol to two 
higher excited levels, n and n — 1, is straightforward. In 
this case, the counterintuitive motional sequence of the 
three PT traps should be performed fulfilling conditions 
Qn,s(do)T > 10 and n„_i ia (do)T ^C 1, which generalizes 
expression (I15p . 

In Fig. [3]we plot curves ft n ,s(do)T = 10 (solid blue) 
and Q n ^i iS (do)T = 1 (dashed red) in the parameter plane 
(do / a, uj x T) for the potential depth s — 6. Perfect fil- 
tering could be performed for each n in the correspond- 
ing grey region transferring the population of state n to 
the right trap without modifying the population of state 
n — 1. From Fig. [3] it is clearly shown that the opti- 
mal minimum distance for the filtering protocol decreases 
with n. This minimum distance, <i™ ln , for the filtering 
protocol involving states n and n—1 can be estimated as 
follows. From Eq. (11) and extending Eq. (15) to excited 
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FIG. 3: (Color online) Curves n„, s (do) T = 10 (solid blue) 
and f2 n _i jS (do)T — 1 (dashed red) in the parameter plane 
(do/a, WiT) for PT potentials with depth s — 6. Filtering for 
the vibrational level n can be achieved in the corresponding 
grey region. 



states, one obtains: 
O (d u 



A n , s e d '™ > 10, 



Hn-l,. « in ) 

and taking the lower limit of Eq. (17) 

10 , 



ln( 



A-r. 



(17) 



(18) 



d™ ln almost perfectly matches the minimum distance do 
of the lower corner of the n-th grey region in Fig. [3] 



be emptied and the protocol will be resumed. After step 
k = N the population of the left trap is expected to be 
in its ground vibrational state. Therefore, the last step 
k = N + 1 will consist in directly measuring the total 
population at the left trap, i.e., P T [N+ 1] = J2 n P n i N ]- 
At the end, the set {k, P T [k]} will be the result of the 
tomography of the initial population distribution at the 
left trap. To evaluate the efficiency of the quantum to- 
mography protocol we define the following fidelity: 



F QT = 1 . 





E 



P^[k = 0}-P T [k = N+l 



(19) 



where P^ [k = 0] is the initial population of state n in the 
left trap, while P T [k = N + 1 — n] is the total population 
measure at the end of each step k. 

Filtering conditions for the step k of the protocol in- 
volving states n and n — 1 read Cl„ tS (do[k])T > 10 and 
f2„_i, s (<io[fc])T <C 1 implying, as shown in Fig. [3j that 
the minimum distance do [k] at each step k should be de- 
creased. Approximated values for do[k] at each step of 
the protocol could be estimated by using expression (|18|) . 
However, we will use, in what follows, accurate values of 
do[k] by numerically integrating Eq. ([7]). 

To illustrate the technique outlined above, let us con- 
sider three coupled identical PT potentials ((2|) with depth 
Vo = —I56hu} x [s = 12] supporting 12 bound energy lev- 
els. Initially, only the left trap is populated with a trun- 
cated thermal distribution among the lowest eight states 
given by: 



IV. 



MULTIPLE-STATE FILTERING 



P^[k = 0] = Ze: 



withn = 0, ..., 7 



(20) 



In this section, we will extend the previously discussed 
adiabatic passage technique to the situation where the 
atomic population is initially distributed among N + 1 
vibrational states of the left PT trap. We will discuss 
first a detailed protocol to perform quantum tomography 
of the atomic population distribution at the left trap and 
later on we will briefly outline a similar approach for 
quantum engineering of Fock states. 



A. Quantum Tomography 

To perform quantum tomography, we will apply the 
filtering protocol sequentially in N steps, i.e., state by 
state, from the most excited [n = N) down to the first 
excited (n = 1) state. At each step k = {1, ...,N}, we 
will transfer to the right trap the population of the corre- 
sponding excited state n = N + 1 — k and keeping the rest 
(from n—1 to 0) in the left trap. After each step k, the to- 
tal population in the right trap, P T [k] = J2 n Af M; w ^ n 
Pn[k] being the population of state n in the right trap, 
will be computed, i.e., measured from the experimental 
point of view. After the measurement, the right trap will 



where Z is the normalization constant and /3 is propor- 
tional to the temperature. 

Figure|4]shows the distribution of population in the left 
(a) and the right (b) traps at each step k of the quan- 
tum tomography protocol performed with the sequence 
of distances do[k] plotted in (c). The initial distribution 
is depicted in step k = of the corresponding figure be- 
ing the truncated thermal distribution given by Eq. ([20]) 
with /3 = 1 in the left trap (a) while the right is empty 
(b). It is clearly seen that at each step k, the population 
of state n = iV+l — fcis transferred to the right trap 
while the population of states from n - 1 to remains 
in the left trap. The mean and the variance of the pop- 
ulation distribution of the vibrational states of the left 
and right traps at each step of the protocol are shown 
in (c) and (d), respectively. Note that the process leads 
to a sequence of Fock states from N to 1 at the right 
trap, since we assume that after each step the right trap 
is emptied. For the left trap, both the mean value of the 
population distribution and its variance decrease at each 
step giving the ground Fock state after the tomography 
process. The fidelity achieved in the case shown in Fig.0] 
is of F Q T = 0.97. 



(a) Left trap 



(b) Right trap 



(a) Left trap 



(b) Right trap 
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FIG. 4: (Color online) Quantum tomography protocol via a 
sequence of adiabatic passage processes in three PT poten- 
tials. Population distribution for the lowest 8 energy levels 
in the left (a) and right (b) traps after each step k of the 
protocol. The minimum distance between the traps at each 
step, do [k], is shown in (c). Mean (m) (closed squares) and 
variance (An-) (closed circles) of the population of each state 
for the left (i — L) (d) and the right (i — R) (e) traps along 
the protocol. The trap-approaching sequence at each step 
of the protocol is given by Eq. (|13|) with u) x T /2 — 40 and 
vo/acj x = 3 x 10 -3 and the corresponding do[k]. 



B. Quantum engineering of Fock states 

The robustness and selectivity of the filtering protocol 
proposed here, allows us to use it for engineering partic- 
ular Fock states. We have already seen in the previous 
subsection that it is possible to generate Fock states at 
specific excited vibrational levels in the right trap after 
each step of the tomography protocol. Moreover, apply- 
ing the quantum tomography protocol from k = 1 up to 
k = N, corresponding to filtering of states from n = N 
to n = 1, one ends up with the Fock state at the ground 
vibrational level n — in the left trap. We have to men- 
tion that the Fock ground state in the left trap could be 
also reached just by one spatial adiabatic passage process 
with minimum distance d a chosen to perform filtering be- 
tween the ground and first excited states as described in 
Section IIII CI Under such conditions, all excited states 
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FIG. 5: (Color online) Ground state filtering by applying 
a single adiabatic transport process. Population distribution 
for the lowest 8 energy levels in the left (a) and right (b) 
traps before and after the single step filtering protocol. The 
minimum distance do is chosen to fulfill the filtering conditions 
for the first excited vibrational state and it coincides with the 
value used in the last step shown in Fig(4] The fidelity of the 
process is above 99%. Other parameters are the same as in 
Fig. SI 



will be transferred from the left to the right trap simulta- 
neously, since the adiabaticity conditions being fulfilled 
for the first excited vibrational state are fulfilled also for 
all excited states above it. We have checked numerically 
such transfer obtaining fidelities of the process above 99% 
(see Fig©. 

Alternatively, by combining each of the adiabatic pas- 
sage steps with a thermalization process it could be pos- 
sible to implement a tunneling assisted evaporative cool- 
ing protocol. In this case, the protocol would consist in a 
sequence of periodically performed adiabatic spatial pas- 
sage processes but without the need to control precisely 
the minimum distance d,Q. In contrast to the standard 
forced evaporative cooling technique [13( developed for 
magnetic traps, and not easily applicable to dipole traps, 
such tunneling assisted (forced) evaporative cooling could 
be performed quasi-continuously in dipole traps without 
the need to open the trapping potential, which could be 
an advantage in coherent control and coherent manipu- 
lation of trapped cold atoms and molecules. 



V. CONCLUSION 

In this paper, we have addressed the filtering of the 
population of specific vibrational states of a Poschl-Teller 
type potential. For this purpose, we have chained the ini- 
tially populated left trap with two empty identical ones 
and we have performed a vibrational state selective spa- 
tial adiabatic passage process from the left trap to the 
outermost right trap. We have derived analytically the 
filtering conditions for the two-state case either involving 
the ground and first excited states as well as two higher 
excited states and we have applied them to the filtering of 
an arbitrary number of vibrational states. By numerical 
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FIG. 6: (Color online) Gaussian potential (black solid curve) 
approximated with harmonic (red dotted curve) and Poschl- 
Teller (blue dashed curve) potentials. The corresponding 
eigenenergies are represented by horizontal black solid (Gaus- 
sian potential), red dotted (harmonic potential) and blue 
dashed (PT potential) lines. No, Nh and Npt give the 
number of vibrational states for the given potential depth for 
the Gaussian, harmonic and Poschl-Teller potentials, respec- 
tively. 



integration of the Schrddingcr equation, we have demon- 
strated that efficiencies of the protocol above 99% can 
be achieved for a wide set of parameter values leading 
to the transfer of the population of all vibrational states 
above a certain one from the left PT trap to the outer- 
most right trap, while the states below it remain at the 
initial left trap. We have also shown that spatial adia- 
batic passage can be used to perform quantum tomogra- 
phy of the initial population distribution of the left trap 
with fidelities above 97% by applying the filtering pro- 
tocol starting from the most excited state to the lowest 
one and sequentially decreasing the minimum distance 
between the traps at each step. Finally, we have also 
briefly discussed the possibility of quantum engineering 
Fock states and of tunneling assisted evaporative cooling. 



monic ones. In this appendix, we will discuss the con- 
venience to use Poschl-Teller potentials instead of har- 
monic ones to approximate Gaussian potentials. 

A single Gaussian potential in one dimension, black 
solid line in Fig. 5, can be written as 

V G {x) = -V Q cxp(-2x 2 /a 2 ), (Al) 

where Vq denotes the potential depth and a — wq, the 
waist of the focused light beam. The harmonic approxi- 
mation of this Gaussian potential, Eq. (Al), reads: 



V H (x) = -V Q (l-2x 2 /a 2 



(A2) 



and it is depicted in Fig. 5 (red dotted curve). Although 
the eigenenergies and eigenstatcs of the harmonic poten- 
tial can be obtained analytically, it is obvious from Fig. 5 
that the harmonic approximation is only accurate for few 
of the lowest vibrational states. Moreover, the harmonic 
potential has to be cut at some energy value in order to 
give a finite number of energy eigenstates. 

To account for the finite number of bound states of 
the Gaussian potentials it is more convenient to use po- 
tentials with known analytical solutions that support a 
finite number of bound states. One of such potentials 
is the Poschl-Teller (square hyperbolic secant) potential 
Q. The approximation for the Gaussian potential (|A1|) 
with a Poschl-Teller potential is given by: 



V PT {x) 



-Vnsech 



V2x/c 



(A3) 



corresponding to the blue dashed curve in Fig. [6j By 
comparing the three potentials in Fig. [6j it is clear that 
the spectrum of the harmonic potential approximation 
(red dotted curve) gives six equidistant vibrational states 
and that its profile fits only close to the bottom of the 
Gaussian potential. On the contrary, the Poschl-Teller 
potential spectrum (blue dashed curve) consists of twelve 
non-equidistant bound vibrational levels and the approx- 
imation to the Gaussian potential shape is much more 
accurate. 
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Appendix A: Approximations for Gaussian 
potentials 

Optical dipole force potentials created by focused light 
beams exhibit a Gaussian spatial profile. These type of 
potentials, due to the lack of analytical expressions for 
their eigenenergies and eigenstates, are very often ap- 
proximated by other types of potentials, typically har- 



Appendix B: Tunneling rate between two 
Poschl-Teller potential traps 

The Gram-Schmidt (GS) orthonormalization proce- 
dure [12J applied to Poschl-Teller gives accurate analyti- 
cal expressions for the tunneling rate of the ground state, 

Considering the single trap eigenstatcs given by 
Eq. (5), the symmetric </>q" s and antisymmetric </>q~ s or- 
thogonal states for the ground, n = 0, vibrational state of 
two coupled identical PT traps of depth Vo = — s (s + 1) 
are given by: 



ti. 



l 

V2 
(±ch~ s 



(±0o,8 (a) + <t>o,s (x - d)) = 

(V2f)+c^(V2^)) 



V2 



Nj,o,s(Bl) 



One could find the energies Eq s = (6$ s \H\ </> \ 
for the eigenstates (|B1[) and finally the tunneling rate 
^o,s{d) = \Eq s (c[) — Eq s (cL)\ for the ground vibrational 
level of two PT traps separated by a distance d: 



and 



^o,s (d) 



2s 2 W J [e 2a D + e~ 2a E] 
A J 1 - W 2 

[ e 2a(s-l) B + e -2a(s-l) C .j w ' 

l-W 2 



. (B2) 



A = A(s,a) = F x (s;s,s;s + l;-e- 2a ,-e 2a ) , 

B = B(s,a) = Fx(s + l;2,2s;s + 2;-e- 2a ,-e 2a ) , 

D = D(s,a) = Fx(s + l;s,s + 2;s + 2;-e- 2a ,-e 2a ) 



with: 



(,± \a± 



; s |^ s ) = 1±W, 

b ,s(x)\<p ,s(x-d)) = W = N 2 0s -A, 

s 

V2sT(.s + l/2) 



Nlo,s 



V^r(s + i)' 



where Fx (a; /3, f3'; 7; z\, z 2 ) is the Appcll hypcrgeometric 
function, which is a generalization for the hypergeomet- 
ric functions 2-F1 {<x,/3;j;z) to two variables zx and z 2 ; 
C (s, a) = B (s, —a) and E (s, a) = D (s, —a); a = d/y/2. 
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